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ABSTRACT 

The hot- Jupiter WASP- 10b was reported by iMacieiewski et al. I d2011albh to show transit tim- 
ing variations (TTV) with an amplitude of ~ 3.5 minutes. These authors proposed that the 
observed TTVs were caused by a 0.1 Mj up perturbing companion with an orbital period 
of ~ 5.23 d, and hence, close to the outer 5:3 mean motion resonance with WASP-lOb. To 
test this scenario, we present eight new transit light curves of WASP- 10b obtained with the 
Faulkes Telescope North and the Liverpool Telescope. The new light curves, together with 22 
previously published ones, were modelled with a Markov-Chain Monte-Carlo transit fitting 
code. 

Transit depth differences reported for WASP-lOb are thought to be due to star spot in- 
duced brightness modulation of the host star. Assuming the star is brighter at the activity 
minimum, we favour a small p lanet ary radius. We find R p = 1.039t n 'rviQ R.iup in agree- 
me nt with iJohnson et all d2009t) and IMacieiewski et all d201 lbl) . IMacieiewski et alj d2011al) 
and lHusnoo et alj d2012l) find no evidence for a significant eccentricity in this system. We 
present consistent system parameters for a circular orbit and refine the orbital ephemeris of 
WASP-lOb. 

Our homogeneously derived transit times do not support the previous claimed TTV sig- 
nal, which was strongly dependent on 2 previously published transits that have been incor- 
rectly normalised. Nevertheless, a linear ephemeris is not a statistically good fit to the transit 
times of WASP-lOb. We show that the observed transit time variations are due to spot occulta- 
tion features or systematics. We discuss and exemplify the effects of occultation spot features 
in the measured transit times and show that despite spot occultation during egress and ingress 
being difficult to distinguish in the transit light curves, they have a significant effect in the 
measured transit times. We conclude that if we account for spot features, the transit times 
of WASP- 10 are consistent with a linear ephemeris with the exception of one transit (epoch 
143) which is a partial transit. Therefore, there is currently no evidence for the existence of 
a companion to WASP-lOb. Our results support the lack of TTVs of hot-Jupiters reported for 
the Kepler sample. 

Key words: stars: planetary systems - stars: individual (WASP-10) -stars: spots -techniques: 
photometric 



1 INTRODUCTION 



* E-mail:susana. barros@oamp.fr 



The discovery of WASP-lO b, a 3 M Jup hot - Jupite r planet in a 3.09 
day orbit, was reported by IChristian et al.l (2009). Its host star is 
a K5V type with T eff = 4675 ± 100 K, [M/H] = 0.03 ± 0.2 
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dChristianet all [2009) and is relatively young, with an age of 
270 ± 80 Myr dMacieiewski et alj|2011al). Using a high signal-to- 
noise transit light curve, Ijohnson et alT ([2009) updated the stellar 
density and re-derived the stellar mass to A/* = 0.75 ± 0.04 Mq 
and radius to R, = 0.698 ± 0.012 Rq. There has been som e 
discussion about the pla netary size, with IChristian et al.1 d2009l) . 
iDittmann et alj d2010h and lKreicova et al] j2010h obtaining a plane- 
tary radius of 1.28±0.09 R^JwimeEotoso^^lJiooi) derived 
a smaller radius of 1.08 ± 0.02 Rj up . iMacieiewski et al.1 J201 lbf) 
argue that the discrepancy in previous results was due star spot in- 
duced brightness variability. In fact, WASP- 10 has been shown to 
have stellar variability due t o the rotation mod ulation of spots with 
a period of 11.91 ± 0.05 d dSmith et alj2009h and semi-amplitude 
up to 10.1 mmag. Star spots reduce the effective stellar disc area 
leading to an underestimation of the stellar radius and a larger tran- 
sit depth. This variation of transit depth with spot cover age explains 
the di fferences in the measured planetary radius (e.g. ICzesla et al.l 
2009). Using a transit taken close to m aximum brightness t hat cor- 
responds to a minimum spot coverage, Macieie wski et al] d201 lbl) 
obtained a maximum planetary radius of 1.03lg Q3 Rj up . 

Interestingly, IMacieiewski et al.1 d201 lal) reported a periodic 
modulation of transit timing variations (TTVs) for WASP- 10 with 
an amplitude of ~ 3.5 ± 1.0 minutes. They suggested a 0.1 Mj up 
planet companion in the outer 5:3 me an motion resonant orbit as 
the cause of the measured TTVs. Later. IMacieiewski et al.ld201 lbl) 
confirmed the TTV solution and showed that only spots located 
close to the stellar limb would affect the transit times but they 
would also result in transit duration variations. Since no transit du- 
ration variations were found in their sample, the authors argued that 
spot occultation features did not significantly affect the measured 
transit times of WASP- 10. 

The possibility of finding additional planets by measuring 
their effect on the centra l transit times of known t ransi ting planet 
system s was proposed bv lHolman & Murravl (2005 ) and Ag ol et al.l 
(2005). They showed that for systems near mean-motion reso- 
nances this method can be sensitive to planets less massive than 
the Earth. This has led many groups to search for TTVs from the 
ground. However, most of the searches resulted in only upper lim- 
its on the mass o f possible companions (e. g. [Miller-Ric ci et al] 
120081 : lBeanll2009l ; iGibson et al.1 120091 , l20ld) . The few initial re- 
ports of TTVs were unconfirmed later. For example, iDfaz et ail 
d2008h found TTVs in OGLE-l l lb which were unconfirmed in 
later studies d Adams et al. l l20ld : lHoveretalj|20l 1). At present, 
only three hot- Jupiter transiting systems s how TTVs that were mea - 
sured from ground-based data: WASP -3 dMacieiewski et al1l2010h. 
WAS P- 10 iMacieiewski et al] 1201 lal) . and WASP-5 iFukui et al] 
2011). Meanwhile, Kepler dat a has allowed the dis covery of many 
TTV systems (e.g. K epler-9 dHolman et al] |20 loh and Kepler- 11 
dLissauer et alj|201lt) ). However, a recent TTV study of the Ke- 
pler sample revealed the lack of TTVs in hot-Jupiters, contrasting 
with the TTVs fou nd for longer-period Jupiters or hot-Neptunes 
dSteffen et all2012l) . The authors suggested that this implied a dif- 
ferent formation and/or evolution history for hot-Jupiters. There- 
fore, it is very important to confirm and understand the hot-Jupiter 
TTV detections that have been reported from ground-based data. 

To confirm and possibly characterise WASP-lOb's companion 
we obtained new high precision transit light curves of WASP- 1 0b. 
Six transit observations of WASP- 1 0b were obtained with the RISE 
fast CCD camera mounted on the Liverpool Telescope and two 
more with the Merope camera on the Faulkes Telescope North. We 
combined the new observations with previously published 22 light 
curves of WASP- 10 to homogeneously derive the transit times. 



The new observations are described in Section 2. In Section 
3, we discuss our transit model. We present the updated parameters 
of the system in Section 4 and in Section 5 discuss our results and 
show possible causes of spurious TTVs. Finally, we present our 
conclusions in Section 6. 



2 OBSERVATIONS 

We present further high precision transit observations of WASP- 
10b taken with Faulkes Telescope North (FTN) and with the Liver- 
pool Telescope (LT). 

Two full transits of WASP-lOb were observed on the 2008-09- 
15 and 2008-09-19 with the Merope camera on the FTN using the 
SDSS-i filter. The Merope camera consists of a e2V, 2048 x 2048 
pixel-CCD and has a field of view of 4.7' x 4.7' . Ap erture photom- 
etry was performed with the DAOPHOT package dStetson|[l987h 
within the IRAF4 environment. The differential photometry was 
performed relative to at least 5 comparison stars that were within 
the FTN field-of-view. 

Anoth er six transits we r e observed with th e fast CCD cam- 
era RISE dSteele et al.ll2008l : IGibson et alj [2008) mounted on the 
2.0m Liverpool Telescope on La Palma, Canary Islands. RISE 
has a wideband filter covering ~ 500-700 nm which corresponds 
approximately to V+R. The pixel scale is 0.54 arcsec/pixel re- 
sulting in a 9.4' x 9.4' field-of-view. Thanks to its frame trans- 
fer CCD, RISE has a deadtime of only 35 ms for exposures 
longer th an 1 second. To d ecrease systematic noise due to poor 
guiding dBarros et al] 1201 ll) the telescope was defocussed to - 
0.6mm and we used an exposure time of 9 seconds. However, 
in one of the observations (2010-11-10) we experimented with a 
higher defocussing of -1.0mm and longer exposure time of 34 sec- 
onds. The RISE dat a were reduced using the ULTRACAM pipeline 
dDhillo n et al. 2007) following the same procedure as for WASP- 
21b dBarros et al] 201 ll) . Details of the observations are given in 
Tabled 

The final new high precision transit light curves for WASP- 
10b are shown in Figure [T] along with the best-fit model described 
in Section [33] We overplot the model residuals and the estimated 
uncertainties, which are discussed in Section |3~!2l 

Previously published WASP- 10b light curves were included 
in our analysis which were kindly made available by the authors. 
Details of these observations are given in Table [2] For complete- 
ness we also show the new best-fit model overploted on the previ- 
ously published transit light curves in Figure [2] whose references 
for these data are given in Table [2] 

Hereafter, we refer to each light curve by their epoch num- 
ber relative to the ephemeris presented in Section 14.21 The epoch 
number is also given in the observation logs and light curve plots. 



3 DATA ANALYSIS 

3.1 Re-analyse of the radial velocity observations 

The discovery paper for WASP- 10b dChristian et al.l2009l) favoured 
an eccentric orbit for the planet e = 0.059^p nnl and lu = 
2.917^0 345). However. IMacieiewski et al] d201 lah argued that ec- 
centricity might have been overestim ated due to stellar v ariability 
and favoured a circular orbit. Later, iHusnoo et al.l d20 1 2h showed 
that the eccentricity detection was dependent on the first two ra- 
dial velocity (RV) SOPHIE observations that could have a dif- 
ferent zero point. Using a method of Bayesian model selection 
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Figure 1. Phase-folded new light curves for WASP- 10. From top to bottom and left to right in chronological order; FTN2008 September 15 and 19, RISE 2008 
October 04, RISE 2010 September 03 and 09, RISE 2010 October 13, RISE 2010 November 10, RISE 2010 December 14. The light curve name including 
the epoch number is printed in each plot. For each light curve, we superimpose the best-fit transit model 2 (solid red line), show the residuals from this model 
at the bottom of each the plot, and give the estimated red noise parametrised by j3 inside parenthesis. For some light curves we also show a transit model that 
includes a spot feature (dashed green line). The light curves and models have been shifted to the TO determined with model 2. The RISE data are binned into 
45 second periods except for the 2010 November 10 which had a longer exposure time. The individual light curves plotted here are available in electronic form 
at CDS. 
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Figure 2. Phase-folded light curves for WASP- 10 from previously published light curves which are described in Table [2] Other comments are the same as 
Figure[T]but red noise values equal to 1 were omitted. 
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Figure 2 - continued 
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Table 1. WASP- 10b new observations log. 



Epoch Date telescope exp time number exposures aperture radius 







sec 




arcsec 




20 


2008-09-15 


FTN 


125 


115 


6.1 


21 


2008-09-19 


FTN 


125 


109 


7.0 


26 


2008-10-04 


RISE 


5 


1580 


8.6 


252 


2010-09-03 


RISE 


9 


1195 


8.1 


254b 


2010-09-09 


RISE 


9 


1233 


8.6 


265 


2010-10-13 


RISE 


9 


1063 


10.3 


274 


2010-11-10 


RISE 


34 


339 


11.4 


285 


2010-12-14 


RISE 


9 


755 


9.2 



Table 2. WASP- 10b previ ous observations log. The reference s are A = IChristian et alj d2009l) . B = I Johnson etail 12009)) . C = iDittmann et all l2010l) . D = 
iKreicova et ai] <2010h , E = lMacieiewski et alj (201 lal) . and F = lMacieiewski et alj l20lTbh , 



Epoch 


Date 


telescope 


exp time 


num exp 


filter 


ref 








sec 




arcsec 




-103 


2007-09-01 


Mercator 


30 


170 


V 


A 


-89 


2007-10-15 


Tenagra 


125 


125 


I 


A 


-88 


2007-10-18 


Tenagra 


125 


91 


I 


A 





2008-07-16 


2.2m Hawaii 


50 


194 


sloan z' 


B 


36 


2008-11-04 


Slovak 




120 


R 


D 


56 


2009-01-05 


0.6m Torun 


50 


254 


R 


E 


123a 


2009-07-31 


0.6m Jena 


70 


144 


R 


E 


123b 


2009-07-31 


0.6m Rozhen 


120 


116 


R 


E 


124 


2009-08-03 


0.6m Rozhen 


90 


122 


R 


E 


132 


2009-08-28 


0.6m Rozhen 


90 


99 


R 


E 


133a 


2009-08-31 


Slovak 




165 


R 


D 


133b 


2009-08-31 


0.6m Jena 


55 


203 


R 


E 


134 


2009-09-03 


0.6m Rozhen 


90 


105 


R 


E 


143 


2009-10-01 


2.0m Rozhen 


15 


548 


V 


E 


144 


2009-10-04 


Slovak 




118 


R 


D 


145 


2009-10-07 


Slovak 




114 


R 


D 


147 


2009-10-14 


Kuiper 


10 


579 


I 


C 


158 


2009-11-17 


Tenagra 


76 


76 


R 


E 


242 


2010-08-03 


2.2m Calar Alto 


40 


165 


R 


F 


243 


2010-08-06 


2.2m Calar Alto 


40 


213 


R 


F 


253 


2010-09-06 


2.2m Calar Alto 


50 


164 


R 


F 


254a 


2010-09-09 


2.2m Calar Alto 


60 


184 


R 


F 



they find the eccentricity detection is not significant and obtained 
e = 0.052 ±0.031. Furthermore, the eccentricity is very correlated 
with the offset between the SOPHIE and FIES data. For all these 
reasons we assume a circular orbit. A better constraint on the ec- 
centricity requires further RV observations. To estimate the stellar 
reflex velocity for a circular orbit we re- analysed the RV observa- 
tions presented in lChristian et al 1 d2009h using the new ephemeris 
given in Section l4~2l 

Using an MCMC routine we fitted 5 parameters of a keple- 
rian model for the host star reflex motion: orbital period, central 
transit time (TO), RV semi-amplitude (K), centre-of-mass velocity 
(7), and an offset between SOPHIE and FIES data. We imposed a 
Gaussian prior to the orbital period and TO using the new ephemeris 
given in Section 1421 and obtained K = 506.4 ± 7.4ms -1 , 7 = 
-11427.9 ± 5.3 ms" 1 , and offset= 112.0 ± 1.2 ms" 1 . This so- 
lution has median rms of 43 m/s which compares with the median 
RV uncertainty of 20 m/s. Hence, the radial velocities have extra 
jitter which could be due to stellar activity. 

iMacieiewski et al] d2011al) argued that the stellar activity had 



a sinusoidal effect (semi-amplitude of 65 m/s) in the measured 
RVs which should be accounted for when deriving the stellar 
reflex velocity due to W ASP- 10b. Following the procedure of 
Macieiews ki et alj ( l2011al) . we computed the Lomb-Scargle peri- 
odogram of the RV residuals after removing the circular model 
given above. This is shown in Figure [3]a s well as the stellar rotation 
period of 11.91 ± 0.05 d determined bv lSmifh et all J2009t) and the 
false alarm probability at 2 a (95.4%) confidence limit. We find no 
significant peak at the stellar rotation period, and hence we cannot 
confirm that the radial velocities are modulated by spot variability. 
Moreover, other host stars that show 1-2% photometric variability 
have smaller activity induced jitter in the measured RVs than the 
claimed 65 m/s for WASP- 10. F or example, HP 18 9733b has a RV 
jitter semi-amplitude of 15 m/s feoisse et alj|2009h a nd CoRoT-7b 
has a RV jitter semi-amplitude smaller than 20 m/s dOueloz et al.1 
2009). Therefore, we conclude that the extra jitter observed in the 
RVs of WASP-10 does not appear to be periodic. Including one 
jitter parameter for each data set in the radial velocity fit we find 
that the FIES data jitter is less than 13 ms" 1 while the SOPHIE 
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Table 3. Estimated red noise values parametrised as j3 calculated according 
to lWinn et aTU2008h 



Epoch beta ind beta simu 



Figure 3. Lomb-Scargle periodogram of the RV residuals after subtracting 
a circular orbit for WASP-lOb. The vertical line shows the rotation period 
of 11.91 and the horizontal line shows the false alarm probability at 2 a 
confidence limit. 



data jitter is 41 ± 17 ms . The latter is much hig her than the 
expect ed jitter for SOPHIE supporting the result of iHusnoo et al.l 
about the possible offset of the two first SOPHIE points rel- 
ative to the subsequent observations. Thus, we experimented dis- 
carding the first two SOPHIE points from the fit which resulted in 
a much lower estimated jitter for the SOPHIE data = 13±18 ms -1 
compatible with the expected values for this instrument. This was 
adopted as our final solution for the RVs with the remaining pa- 
rameters being: FIES jitter 8 ± 13 ms" 1 , K — 540 ± 11 ms' 1 , 
7 = -11461 ± 13 ms" 1 , and offset= 84 ± 18 ms" 1 . Our final 
solution for t he stellar reflex veloci ty is in agreement with the value 
presented by Christian et aL (2009). The expected RV signal of the 
possible companion o f 14 m/s is below the media n RV uncertainty, 
as also pointed out by Macieiews ki et al.l feOllal) . 

3.2 Photometric errors 

An accurate estimate of the photometric errors is important to ob- 
tain reliable system parameters. We begin by scaling the errors of 
each light curve so that the Xred °f the best fitting model is 1.0. 
This is especially important in our case because we are combining 
data from different instruments which were reduced with different 
pipelines. 

Exoplanet transit observations are often affected by correlated 
red noise which must be taken into account in error analyses. To 
estimat e the time-c orrelated noise factor f3 we followed the proce- 
dure of I Winn et al.l J2008h . This consists of comparing the standard 
deviation of the residuals of the best fit solution for different time 
bins. For each light curve, we rescaled the errors by the estimated 
P which is given in Table [3] This procedure does not distinguish 
between systematic noise and real stellar variability. All departures 
from the transit shape are considered as noise. For example, when 
we fit simultaneously all light curves (model 2), variations of the 
transit shape will also be considered as noise. Therefore, the values 
of /3 for the simultaneous fit are higher than for the individual fits. 
If the shape variations are real this procedure will overestimate the 
errors. However, if they are indeed caused by systematic noise it 
will allow a more accurate error estimation. 



3.3 Transit model and fitting procedure 

To determine the planetary and orbital parameters, we fitted the 
light curves of WASP-lOb with a transit model coupled with an 
Mark ov Chain Monte Ca rlo (MCMC) routine, a procedure simi- 
lar to iBarros etaL (201 ]]) for WASP-21b. In our previous transit 



-103 
-89 

-88 


20 
21 

26 

36 

56 

123a 

123b 

124 

132 

133a 

133b 

134 

143 

144 

145 

147 

158 

242 

243 

252 

253 

254a 

254b 

265 

274 

285 



1.0 
1.0 
1.0 
1.0 
1.0 
1.5 
1.5 
1.4 
1.0 
1.1 
1.1 
1.1 
1.0 
1.4 
1.7 
1.3 
1.5 
1.4 
1.9 
2.4 
1.0 
1.5 
1.8 
2.0 
1.7 
1.0 
3.0 
2.7 
1.8 
1.3 



1.0 
1.0 
1.0 
1.0 
1.3 
1.7 
2.0 
1.4 
1.0 
1.1 
1.1 
1.0 
1.0 
1.7 
1.8 
1.3 
1.8 
1.6 
2.3 
4.3 
1.0 
1.8 
2.3 
2.0 
2.0 
1.2 
3.0 
3.0 
2.2 
1.5 



fitting analyses we fit all the light curves simultaneously. How- 
ever, here we also present the individual light curve fits, mostly 
for comparison with previous result s and call the latter mo del 1 
and the former model 2. We used the iMandel & Agoll 00021) tran- 
sit model parametrised by the normalised separation of the planet, 
a/R*, ratio of planet radius to star radius, R p /R t , orbital incli- 
nation, i, and the transit epoch, To, of each light curve. A circu- 
lar orbit was adopted as discussed above. For each light curve, we 
also include two extra parameters (F ou t and T gra d) to account for 
a linear normalisation. We used the quadratic limb darkening co- 
efficients (LDCs) deduced from the tables of iHowarthl d201l|) for 



:ff 



4675 K, log g= 4.4 and [M/H] = 0.03 ([christian et al. 



2009). Most of the light curves are of insufficient quality to allow 
for the fitting of the LDCs. To be consistent in the treatment of all 
the light curves we fixed the limb darkening coefficients to the the- 
oretical values. Although this may lead to a slight underestimation 
of the errors, it was preferred to avoid spurious correlations. The 
values of the LDCs used for each filter are given in Table [4] 

The MCMC chains were constructed following the 
Metropolis-Hastings algorithm. At each step of the chain 
each parameter is perturbed by a jump function, which is a 
random number sampled from a Gaussian distribution with mean 
zero and standard deviation equal to the respective parameter 
uncertainty. The new parameter set is accepted with probability: 

P — rain (l, exp ^ j) where A\ 2 is the difference in the 

X 2 of subsequent parameters sets. The jump functions are scaled 
by a common factor in order to ensure that the percentage of 
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Table 4. Limb darkening coefficients used for each filter. 



filter 


linear coef a 


quadratic coef b 


sloan z' 


0.391 


0.210 


I 


0.421 


0.192 


R 


0.602 


0.133 


rise 


0.612 


0.146 


V 


0.763 


0.036 



accepted steps lies between 20% and 30 %. Further d etails about 
the fitting procedure can be found in lBarros et al.l (1201 II) . For both 
models the initial 20% of each MCMC chain that corresponded to 
the burn-in phase were discarded and the remaining parts merged 
into a master chain. The best fit parameter was estimated as the 
median of its probability distribution and the 1 a limits as the value 
at which the integral of the distribution equals 0.341 from both 
sides of the median. To test convergence, the Gelman & Rubin 
( 1992) statistic was computed for each fitted parameter and was 
found to be less than 0.9% from unity for all the parameters in 
both models. 



3.3.1 Model 1 

The existence of an additional planet or moon affects the tran- 
sit timings of the transiting planet, but can also affect the or- 
bital parameters of the system. Hence, in previous analyses of 
the transit times of WA SP-10, the authors op t ed for fitting each 
light c urve individually (Dittma nn et al.l [2010l ; iMacieiewski et al.l 
1201 1 al ibi). However, in some cases they were forced to fix the in- 
clination of the system due to the poor quality of individual light 
curves. Here, we follow the same procedure primarily for compari- 
son with previous methods. The results will also test possible vari- 
ations of the light curve shape. Thus, we fitted each light curve in- 
dividually using the MCMC method explained above with 500 000 
points per chain. Five chains were combined to obtain the final re- 
sult. Instead of fixing the inclination we choose to impose a prior 
of the form: 



where we assume in = 88.66, <n = 0.12, taken from the simul- 
taneous fit (model 2). This allows to account for the uncertainty of 
the inclination in the fit. We refer to this as model 1 . 



3.3.2 Model 2 

For model 2 we assume that a possible additional planet has no sig- 
nificant effect on the transit shape during the time span of the obser- 
vations. Therefore, we performed a simultaneous MCMC fit to all 
30 light curves, fitting globally a/ R. t , R p /R,, and i and individu- 
ally TO, F ou t and T gra d- We remind the reader that deviations from 
a constant shape are considered to be only due to systematic noise 
or intrinsic stellar variability and hence are treated as noise. Specif- 
ically, deviations of limb darkening parameters from the tabulated 
values are also considered as noise. However, since the brightness 
of WASP-10 has been shown to vary s ignificantly due to star spots 
dSmifh et af]|2009l ; IMacieiewski et alj|20lTah a factor df was in- 
cluded to account for out-of-transit flux variability between each 
light curve. This factor has b een set arbitrari l y to z ero for our best 
quality light curve (Epoch 0, Ijohnson et al.l d2009h ) which is also 



the one taken in the reddest filter, and hence that should be less af- 
fected by stellar variability. Therefore, a total of 30 x 4 + 3 = 123 
parameters were fitted simultaneously. We computed seven MCMC 
chains each of 1 500 000 points and different initial parameters. 



4 RESULTS 
4.1 Model 1 

The individual light curve fitted parameters are given in Table [5] 
Transit times were converted to Barycentric Dy namical Time 
(TDB ) using the IDL codes kindly made available bv lEastman et al.l 
(2010). The weighted mean of the fitted parameters is: a/R, — 
11.910 ± 0.024, R p /R, = 0.15869 ± 0.00019, and i = 88.67 
degrees and the dispersion of inclination is equal to 0.13 degrees 
as expected from the prior. 

Comparing the values of a/R* for all light curves we con- 
clude that individual deviations from the mean value of the nor- 
malised separation of the planet are smaller than 2a except for 2 
transits. These are epoch 20 which was observ ed in 2008-09-15 
with FTN and epoch 147 dPittmann etailbOlCl) . Both have a/R, 
higher than the mean a/ R, by more than 3a. Both of these transits 
show signs of correlated noise either due to systematics or stellar 
variability and will be discussed in section 5.2. 

As expected, due to spot induced brightness variability of the 
host star, the parameter R p / R* shows a wide spread. Two of the 
light curves (243, 254a) show a significantly higher (3 a) value of 
Rp/R,, while three light curves (0,133,158) show a significantly 
lower value. Therefore, we confirm transit depth variations reported 
by other authors and the need to include the extra df parameter in 
the global fit. For completeness we also tested a model with free 
inclination and concluded that the inclination is constant within 
the uncertainties. Hence, we conclude that there are no significant 
shape variations apart from the transit depth, except for transits 
epoch 20 and 147. 

Besides testing variations of transit shape, model 1 is also 
useful for direct comparison with previous transit times. This al- 
lows to better isolate the causes for the differences in the estimated 
transit times. The first three transits (-103,-89,-88) were previously 
combined together with the WASP data in a single ephemeris, and 
hence, they cannot be directly compared. Regarding the other previ- 
ously published transit times, our results are within la of previous 
results except for epochs 36 and 147 for which the difference is re- 
spectively — 2.0±0.94 and 0.67±0.47 minutes. The disagreement 
can be explained by failure to account for the baseline function in 
the previous analysis of these transits. In our model a linear trend is 
fitted simultaneously with the transit model. When we do not fit a 
trend we recover the previous published results. TO correlates with 
the baseline function because the latter affects the symmetry of the 
light curve. Therefore, including the baseline function is crucial to 
obtain reliable TOs 

In general, we obtain slightly larger errors (~ 1.6 x) than pre- 
vious studies. For each case this is due to one or more of the fol- 
lowing differences from the previous results: inclusion of a base- 
line function in the fit, inclusion of the red noise factor, and/or the 
treatment of the inclination. We argue that fixing or artificially con- 
straining the inclination will cause underestimation of the transit 
time errors. 

Despite the difference from previous results being small, the 
new transit times are in much better agreement with a linear 
ephemeris. The \ 2 °f a linear fit to the transit times is 91.95 for 
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Figure 4. Lomb-Scargle periodogram of residual transit times from model 
1 after subtracting a linear ephemeris. The vertical dash-dotted line shows 
the rotation period of 11.91. The dashe d lines show the TTV frequencies 
proposed by Macieiewsk i et al.l {201 lal) and the horizontal line shows the 
false alarm probability at 1 and 2 <r confidence limit. 



28 degre es of freedom, i.e. Y rprl = 3 .28 which compares with 13.8 
found by Macieiews ki et al. I feOllal) . However, the p-value is still 
low 9.777 x 10 -9 and hence the null hypothesis can be discarded 
at 5.7 a. Four of the new transit times deviate more than 3 a from 
a linear ephemeris. These are epochs 20, 21, 143, and 147. Both 
epoch 21 and epoch 147 show trends in the light curve that can be 

affecting the transit times. 

To test the model for the TTVs proposed by Macie iewski et alj 
(2011a) we fitted the sin usoidal model given by equation 1 of 
IMacieiewski et alJd201 1 <J) for frequencies: f{ tv = 0.183 cycl P' 1 



and fi 



0.175 cycl P , where P is the orbital period of 



WASP- 10b. The preferred solution presented by Macieiews ki et al.l 
d201 lah for frequency /**" had a p-value equal to 3.4 x 10 -4 which 
suggests that the model could be rejected at the 3.6<r level. How- 
ever, we remind the reader that this sinusoidal model is only used 
to estimate the most probable TTV frequency to be used later in 
the 3-body simulations. Therefore, the p-value test cannot be di- 
rectly applied. Nevertheless, we found that the sinusoidal model 
fit to the new transit times is worse, as it results in a much lower 
p-value 4.1 xl0~ 7 for both frequencies f\ tv and /"". Most impor- 
tantly, for the new results we obtain 3.5 times smaller TTV semi- 
amplitudes: A[ tv = 0.48 ± 0.18 and Af v = 0.53 ± 0.14. 

In Figure [4] we show the Lomb-Scargle periodogram of 
the residual transit times of model 1 after subtracting a lin- 
ear ephemeris. Notice th at the TTV frequencies proposed by 
IMacieiewski et al] J2011ah have disappeared which agrees with the 
much lower amplitudes found for the sinusoidal model fit. 



4.2 Model 2 

The individual transit times and the out-of-transit flux variability 
derived from a simultaneous fit of the 30 light curves of WASP- 10 
are given in Table [6] 

We conclude that three light curves 242, 243, 254a 
( IMacieiewski et al1l2011bh have a significantly lower df (at 3<r ) 
than the chosen zero point. In our notation a negative df means the 
star is fainter which corresponds to a larger transit depth. Some 
light curves have a higher df, but this is non-significant. There- 



Figure 5. Comparison between our results and previously published ones. 
We show the time difference in the top panel and the error ratio in the bottom 
panel. Model 1 comparison is shown as black stars and model 2 as red 
circles. 



fore, we confirm that our chosen zero point is valid. However, this 
zero point might not correspond to the real maximum brightness of 
WASP- 10, and hence the planet to star ratio could still be overesti- 
mated. 

In Figure[5]we plot the difference between the our transit times 
derived with model 1 (black stars) and model 2 (red dots) and previ- 
ously published transit times. We also show the ratio of transit time 
errors. The highest difference in times is found for transits 36, 132 
and 147. For epochs 36 and 147 this is due to the baseline function 
not been accounted for in previous fits, as mentioned before. For 
epoch 132 the difference is within the errors and it is due to epoch 
132 being a partial transit and hence very sensitive to the details of 
the model fit. The transit times derived from model 1 and model 
2 are consistent within 1 a. The biggest differences in the results 
from the two models are in partial transits. In model 2, the global 
shape of the transit defined by the complete transits constrains the 
shape of partial transits, this allows a better constrain on the base- 
line function that correlates with the transit times. 

From the bottom panel of Figure [5] it is also clear that our 
estimated errors are higher than previous ones. This is due to the 
inclusion of red noise in our analysis. For model 2 the uncertainties 
of TOs of partial transits are smaller than for model 1 because the 
shape of the transit in model 2 is almost fixed by the other light 
curves. In contrast, some TO uncertainties of complete transits are 
higher because differences from the common shape are considered 
red-noise which is accounted for as explained in Section l3~2l This 
is the key difference between the models and implies that in this 
case model 2 will produce much more reliable results, as will be 
discussed later. 

We compute a linear ephemeris using the new estimated transit 
times for model 2: 



T t (TDB) = T(0) + EP, 



(2) 



where P = 3.09272932 ± 0.00000032 and T = 
2454664.038090 ± 0.000048. Time residuals (or TTVs) from the 
linear ephemeris are given in Tableland shown in Figure [6] The 
RISE data are show as stars, the new FTN data as squares, and the 
new results based on previously published light curves as circles. 

The x 2 °f the linear ephemeris is 86 for 28 degrees of freedom 
and the p-value is 5.7 x 10 -8 so we conclude that the linear fit 
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Table 5. WASP-10 MCMC fitted parameters for model 1 (individual fit). Times are in TDB. 
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12 
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12 
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1 1 
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n 
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70 


1 2 
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n 
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12 
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21 


12 
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12 
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96 


1 1 
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11 
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56 
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Figure 6. Model 2 estimated transit residuals of WASP-lOb after subtracting a linear ephemeris. The transit times were derived from a simultaneous fit of all 
the light curves. Previously published fitted light curves are shown as circles. The new RISE data are shown as stars while the new FTN data as squares. 



is not a a good fit at 5.4a confidence level. However, only epoch 
143 deviates significantly (at 5.8 a) from the linear ephemeris. This 
transit is the only one obtained with the 2.0m Rozhen telescope and 
is a partial transit. If we discard epoch 143 the \ 2 decreases to 53 



for 27 degrees of freedom and the p-value increases to 0.00176 and 
hence we can only reject the null hypotheses at the 3.1a level. 

For these new transit times deri ved with model 2 we also 
tested the sinusoidal model proposed by Maciejewski et a 
for the two proposed frequencies as done in the previous section for 
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Table 6. WASP- 10 MCMC fitted parameters for model 2 (simultaneous fit). We also give the time residuals from a linear ephemeris in seconds. 



Number 


df 




TO 
days 




Time residuals 
sec 


-103 


0.047 ± 0.065 


4345 


48665 


± 


0.00038 


-28.4 ± 32.6 


-89 


-0.067 ± 0.030 


4388 


78596 


± 


0.00041 


66.7 ± 35.7 


-88 


-0.040 ± 0.057 


4391 


87917 


± 


0.00054 


108.0 ± 46.4 





0.000 ± 0.000 


4664 


03808 


± 


0.00005 


-2.4 ± 4.6 


20 


0.030 ± 0.032 


4725 


89179 


± 


0.00038 


-77.9 ± 32.4 


21 


0.007 ± 0.020 


4728 


98421 


± 


0.00043 


-104.6 ± 37.3 


26 


0.065 ± 0.024 


4744 
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± 


0.0002 1 


13.9 ± 18.3 


36 


-0.066 ± 0.053 


4775 


37696 


± 


0.00054 


51.0 ± 46.7 


56 


-0.029 ± 0.021 


4837 
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± 


0.00043 


5.8 ± 37.4 


123a 


-0.002 ± 0.028 


5044 


44407 


± 


0.0003 1 


20.8 ± 26.9 
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± 


0.00046 


A A 1 I T f\ "I 
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± 
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0.034 ± 0.044 
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± 
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133a 
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5075 


37113 


± 


0.00073 


0.9 ± 62.7 


133b 


-0.022 ± 0.025 


5075 


37133 


± 


0.00043 


18.2 ± 37.4 


134 


-0.028 ± 0.036 


5078 
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± 


0.00062 


90.3 ± 53.5 
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0.023 ± 0.020 


5106 


29655 


± 


0.00032 


-161.2 ± 27.9 


144 


-0.063 ± 0.053 


5109 


39103 


± 


0.00095 


-10.2 ± 82.0 


145 


-0.110 ±0.069 


5112 


48544 


± 


0.00110 


135.1 ±95.4 


147 


-0.026 ± 0.029 


5118 


66795 


± 


0.00056 


-H9.4±48.8 


158 


0.057 ± 0.020 


5152 


68916 


± 


0.00044 


-16.6±37.7 


242 


-0.043 ± 0.014 


5412 


47858 


± 


0.00028 


-4.0 ± 24.2 


243 


-0.053 ± 0.013 


5415 


57163 


± 


0.00026 


23.4 ± 22.4 


252 


0.009 ± 0.021 


5443 


40631 


± 


0.00027 


33.5 ± 22.9 


253 


-0.018 ±0.012 


5446 


49884 


± 


0.00021 


I6.l ± 18.4 


254b 


-0.010 ± 0.021 


5449 


59123 


± 


0.00031 


-12.8±27.0 


254a 


-0.040 ± 0.006 


5449 


59139 


± 


0.00009 


0.7 ± 8.1 


265 


-0.025 ± 0.026 


5483 


61125 


± 


0.00044 


-13.7 ± 38.3 


274 


-0.025 ±0.011 


5511 


44562 


± 


0.00025 


-29.9 ±21.8 


285 


0.037 ± 0.025 


5545 


46534 


± 


0.00029 


-56.4 ± 25.2 



model 1 and found similar results. We find Af v = 0.52 ± 0. 17 and 

A ttv _ q 45 ± q 14 whi , e both p . values Me . less men I g x 1Q -6 

and hence the model could be rejected at the 4.8a level. We reit- 
erate that the expected TTV variations are not strictly periodic and 
the p- value can not be directly applied. However, clearly the new re- 
sults do not support the TTV solution claimed bv lMacieiewski et all 
d2011ah . We analyse a perturber model for the our derived transit 
times in Section [5~Tl 

The periodogram of the time residuals relative to a linear 
ephemeris are given in Figure UJ Comparing with the model 1 re- 
sults, it is clear that the peaks around the star's rotation period have 
decreased, implying that model 2 results are less affected by stellar 
variability. 



4.3 System parameters 

Although the main objective of this study is not deriving updated 
parameters for WASP- 10, for consistency with the simulations pre- 
sented in Section 15.11 we derive the system parameters from the 
geometric fit ted paramete r s, give n in the top of Table [7J The same 
procedure as lBarros et al.ld2012l) is followed. 

We estimate the stellar mass from models by interp olating the 
Yonsei-Yale stellar evolution tracks bv lDemaraue et al.l (120041) us- 
ing the effective temperature T e // = 4675 ± 100 K and m etal- 
licity [M/H] = 0.03 ± 0.2 dex from IChristian et alj J2009h and 
our new derived stellar density p, — 2.359ln n i^/Op). We obtain a 
sightly larger stellar density than derived by Ijohnson et alj d2009h 




0.2 

Frequency (cyc/Pb) 



Figure 7. Lomb-Scargle periodogram of model 2 residual transit times after 
subtracting a linear ephemeris. The vertical dash-dotted line shows the ro- 
tat ion period of 11.91 the das hed lines show the TTV frequencies proposed 
by Macieiewski et al. (201 la) and the horizontal lines show the false alarm 
probability at 1 and 2 a confidence limit. 



( p» = 2.20 ± O.O6P0) which results in a slightly lower mass 
0.73 ± 0.1 Mo, but consistent with previous results. 

Combining the stellar mass, stellar reflex velocity estimated 
in Section [3~T1 and the geometric parameters derived with model 
2, we derive the system parameters which are given on Ta- 
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Table 7. WASP-10 system parameters derived from model 2. The first set 
of parameters are derived directly from the transit light curves. 



Parameter 



Value 



Normalised separation a/R* 
Planet/star radius ratio R p /R* 
Transit duration Ti_4 [hours] 
Transit duration Tti [hours] 
Orbital inclination / [degrees] 
Impact parameter b [R t ] 
Stellar density p* [p©] 

Orbital semimajor axis a [AU] 
Stellar mass M, [M Q ] 
Stellar radius R* [R ] 

Planet mass M p [M Jup ] 
Planet radius R p [Rj up ] 
Planet density p p [pj] 



11.895 ±0.083 

n 1 00036 
U.lO(OO_ 00039 

z.zouo_ 0051 

i n 110 0.0020 
l.»liu_ 0021 



88.66 ±0.12 

7+0.021 
-0.025 
3+0.053 
-0.047 



0.277] 
2.359"* 



0.0375 ± 0.0017 
0.73 ±0.1 



+0.028 



0.678 



3.14 ±0.27 



1.039 
2.81 



+0.043 
0.049 
+0.44 
0.33 



ble [7] The parameters were deri ved using the equations of 
ISeager & Mallen-Ornelasl d2003h and lKippingfeoid) . Specifically, 
the Ti_4 refers to the total transit duration and Tti is the dura- 
tion of the sky-projected centre of the pl anet over l apping the stellar 
surface and is defined by equation 15 of Kip ping] J20K 



5 DISCUSSION 

5.1 Additional perturber 

From the analysis of the TTVs, Macieiews ki et alj d201 lah postu- 
lated the presence of a perturber with a mass of ~ 0.1 A/j up on 
a circular orbit with a period of ~ 5.28 days. This orbital period 
places the system close a 5:3 mean-motion resonance. The exis- 
tence of the planetary companion was motivated by a large reduced 
chi-square Xrod = 13.8 associated to a linear fit of the mid-transit 
times. 

In the present study, additional transit observations and new 
estimations of the error bars lead to Xrcd = 3.1 for a lin- 
ear ephemeris. This valu e is lower than the one obtained by 
Macieiews ki et alj d201 lah . but remains significantly larger than 
one, suggesting the exis tence of an additional signal. As in 
Macieiews ki et alj ^201 lah . we make the hypothesis that this signal 
is produced by a planetary companion and try to constrain its mass 
and orbital properties. For that purpose, we performed 24 x 10 nu- 
merical 3-body simulations with 400 values of the companion or- 
bital period P c ranging uniformly between 1.5 Pi, to 10 Pb, 60 val- 
ues of the initial longitude A c (at the time of the first transit Epoch=- 
103) uniformly distributed over 360 degrees, and 100 masses M c 
between 0.01 to 5 Jupiter masses with logarithmic steps. In all sim- 
ulations, the initial conditions correspond to circular coplanar or- 
bits. The best fit to the timing residuals is very isolated in the pa- 
rameter space and gives P c — 19.3942 days, Al c — 0.2457 Mj up , 
and A c = 48 deg with Xrcd = 1-77. The TTV signal produced 
by this set of parameters is superimposed to the timing residuals in 
Fig.! 

Since our best s olution is completely different from that 
of iMacieiewski et al] J201 lah . we analyse its robustness to the 
mid-transit time uncertainties. We generate 100 sets of artifi- 
cial timing residuals with the same 27 Epochs as the observa- 
tions. The value of each point is then taken from a normal dis- 




Figure 8. Comparison between numerical TTVs (solid blue curve) pro- 
duced by an hypothetical perturber and the timing residuals (red dots with 
error bars). 
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Figure 9. Possible parameters of an hypothetical planetary companion 
based on the timing residuals. Each point is the best fit to a mock data set. 



tribution whose mean corresponds to our best solution and width 
equals the uncertainty of that Epoch. To reduce the cpu time, 
we forgo the 3-body simulations in favour of the perturbation 
method ptmeth l lNesvornv & M orbidelli 2008). This is a pow- 
erful analytical method about 10 4 faster than numerical simula- 
tions w hich provides accurate results o utside of mean-motion reso- 
nances dNesvornv & Morbidelli.20081) . It is also the technique used 
by IMacieiewski et alj d20 1 lah to obtain a good estimate of the lo- 
cation of the best parameters before performing numerical simula- 
tions. The best solutions derived for each of the 100 artificial data 
sets are plotted in Fig. [9] Each point correspond to a different data 
set. The most remarkable result is that these best solutions, with re- 
duced chi-squar^]] of the order of or even lower than 1, are widely 
spread in the mass/period diagram with P c between 1.5 and 12 Pt, 
and M c between 0.01 and 5 Mj up . It is thus not currently possible 
to infer any set of orbital parameters from the analysis of the timing 



1 The values of the reduced chi-square obtained with ptmeth are statisti- 
cally lower than 1.77 found with the 3-body simulations. This is expected 
since the model of the perturbation method contains two additional param- 
eters (eccentricity and longitude of periastron) in comparison to the slowest 
Af-body integration method. Moreover, the mock data sets fitted by ptmeth 
are not affected by stellar activity. 
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residuals of WASP-10. This would require lower error bars and/or 
a larger number of observations. 

It can be noted that the most massive companions at a given 
orbital period, in Fig. [9] are those with the lowest eccentricity. This 
is expected since eccentricities increase the orbital interactions and 
hence the TTVs. This property has already been used on different 
systems to set an upper boundary on the mass of a pertur ber as 
a function of its semi-major axis (e.g. JSteffen & Agol 2005) . For 
instance, in the case of WASP-10, from the A^-body simulations we 
obtain a 3 sigma upper mass equal to ~ 0.1 Mj up at the 5.28 day 
period claimed by Macieie wski et al.| d201 lal). This is not stri ngent 
enough to rule out the solution of IMacieiewski et a and 
neither does it confirm it since it lies at the 3 sigma limit and not at 
the best solution. Moreover, it should be stressed that any perturber 
with mass beyond ~ 1 Mj up would have been detected by radial 
velocities. They can thus be discarded. 

This analysis illustrates the difficulties of interpreting the per- 
turbations produced by a n on-transiting com panion in very low sig- 
nal to noise data (see, e.g- lBoue et alj20~12l) . The detection and the 
characterisation of non-transiting planets do require strong TTV 
signa l s as in the cases of Kepler- 19 and KOT872 f Ballard et al.l 
l201ll : lNesvornv et alj20l3) . 



5.2 The effect of spots in measured TTVs 

Spot occultation f eatures have been shown to affect mea- 
sured transit times jMiller-Ricci et al.l 120081: lAlonso et all 120091 ; 
ISanchis-OiedaetalJl201lh . The transit light curves of WASP- 10b 
show features that might be due to the oc cultation of stellar spots 
by the planet. IMacieiewski et alj 1201 lbh shows that the features 
in epochs 242, 243, 253 are compatible with spots along the tran- 
sit cord with an area of approximately 25% -37% that of planetary 
disc. However, the same shapes fe atures can also be caused by sys- 
tematic noise tearros et al.ll201 lh . In this section, we will test the 
effect of spot features in the light curve shape and consequently 
in the fitted parameters, particularly in the measured transit times. 
Since we are analysing data from different instruments, filters, and 
signal-to-noise ratio we will not attempt to test if the "spot features" 
are due to real spots occultations or systematic noise. Instead, we 
assume that independently of its causes, a "spot feature" will have 
similar effect on the shape of the light curve. 

We modelle d the spot feature s as i nverted transits, a procedur e 
similar to lKundurthv et al. | fcOllh and IMacieiewski et all feOllbl) . 
We used the transit model described above, but for the spot features 
we assume no limb darkening and zero planet-spot impact param- 
eter. Values for a and R p were fixed to those of Table [7] and an 
extra factor was included to account for the spot-to-star contrast. 
This is a very simplified spot model, but reproduces the observed 
features quite well. To obtain typical values for the spot parameters 
for our simulations we used the spot seen in light curve epoch 265 
as a template. According to our tests of light curve systematics this 
appears to be a real spot. 

To simulate the effects of a spot in a transit light curve we 
added the model spot to the model light curve of WASP-lOb at 
different in-transit phases. The resulting light curve was then fit- 
ted with a transit model not including a spot to test how the spot 
affects the fitted parameters. Errors were not included in the sim- 
ulated light curves. We also do not ac count for geometrical fore- 
shortening (San chis-Oieda et all 1201 lh . i.e. the same spot ampli- 
tude was used for each phase. The derived parameters as a function 
of central spot time relative to the mid-transit time are shown in 
Figure [TO] where evidently a spot similar to the one in light curve 




Time (hours) 

Figure 10. The figure shows the effects of a spot in the fitted light curve 
parameters as a function of mid spot time relative to mid transit time. The 
vertical lines show the position of ingress and egress. 



epoch 265 is able to create fake transit time variations up to 1.2 
minutes. Clearly, the spots will also significantly affect the derived 
shape parameters: a/R t , R p /R t , and the inclination. For the best 
signal-to-noise light curves, the spot effect will be more than 3 a of 
the errors quoted in Table [5] 

Noteworthily, the effect of the spot is higher when the spot is 
in the limb of the star, although spots at longitudes up to 30 degrees 
from the limb will significantly affect the measured transit times. 
Furthermore, spot features during ingress and egress are harder to 
detect in the light curve both analytically and by eye. This is illus- 
trated in Figure QT] which shows an example of a light curve where 
the planet occults a spot during egress. The combined light curve 
is well fitted by a transit model with a shorter duration and earlier 
egress time which results in an earlier TO. Moreover, in individual 
fits this light curve would appear to have low red noise. 



5.3 Spot modelling 

As seen in the previous section, the light curves that show higher 
TTVs due to spots also show a larger a/ R, (or shorter duration). 
Hence, a high value a/i?» is an indication that the light curve could 
be affected by spots. Moreover, in this case the red noise might be 
underestimated for individual fits. To check if the TTVs presented 
in Table|6]are affected by spots we tested the light curves that have 
a high red noise and also the ones with high a/ R t . 

From the results in Table|5] as mentioned before, we conclude 
that two light curves have higher a/R, than the mean value. These 
are epochs 20 and 147. Table [3] shows that five light curves have 
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Figure 11. We show a transit light curve with an occulted spot during egress 
as a dashed line. It results from a combination of a transit light curve model 
shown as a solid line and a spot feature shown as a dotted line. Although 
there is no evidence of a spot in the combined light curve, it has a shorter 
duration and a measurable TTV=— 73 seconds. 

red noise higher than 2 for model 2, namely epochs 143, 145, 147, 
254b, and 265. These 6 light curves were fitted with a transit model 
including a spot model as described in the previous section. The 
impact parameter was fixed to and we assume no limb darken- 
ing, the contrast, phase, and duration of the spot were fitted. When 
we include a spot feature we find variations on the fitted transit 
times up to 2.2 minutes (for epoch 20). We extended the analysis 
to all the light curves, but found TO variations only in another light 
curve; epoch 21. The difference between the derived TO when in- 
cluding a spot model and the TO derived with model 2 are given 
in Table [8] as well as the corresponding red noise for the one spot 
model parametrised by f3. For curiosity, we also give the spot fitted 
parameters in Table [8] however the fits are highly degenerate and 
dependent on the starting parameters. We confirm that for all the 
light curves, the time correlated features can be well described by a 
spot, given that the estimate of the red noise significantly decreases. 
However, in some cases such as epoch 147 residual red noise re- 
mains in the light curve after accounting for one spot. Although the 
fits are highly degenerate, we confirm that only spots close to the 
limb of the star have a significant effect on the measured transit 
times. Although the spot induced TTV can be large, the new transit 
times agree with model 2 results. However, this would not be the 
case if we had not included red noise. Furthermore, the existence of 
a spot explains the significant TTVs found for model 1 which does 
not penalise shape variations. Hence, in this case model 1 underes- 
timates the errors in the mid-transit times. 

In Figure [12] we show the timing residuals for the one spot 
model. The x 2 °f a linear ephemeris is 76 for 28 degrees of free- 
dom and the p-value is 1.6 x 10~ 6 which would lead to a rejection 
of the model. Epoch 143 is still inconsistent (5.8 rr) with a linear 
ephemeris since this discrepancy cannot be explained with a spot 
model. The ingress of this transit is clearly earlier than predicted 
by the linear ephemeris which cannot be caused by a spot. If we 
exclude epoch 143 the \ 2 °f me linear ephemeris reduces to 43 
for 27 degrees of freedom (x?ed= 2.56) and the p-value is 0.028 
and hence a linear ephemeris is no longer rejected. Therefore, we 
conclude that if we account for "spot features" the transit times of 
WASP- 10 are consistent with a linear ephemeris. 

Figure[T3]shows the Lomb-Scargle periodogram of the residu- 



Table 8. Variations of the fitted mid-transit times for the one spot model. 
We also show the estimated red noise of the resulting residuals and the spot 
parameters. 
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Figure 13. Lomb-Scargle periodogram of residual transit times from a spot 
model after subtracting a linear ephemeris. The vertical dash-dotted line 
shows the rotation period of 1 1.91 the dashed lines show the TTV frequen- 
cies proposed by Maciejewski et al. 1 201 la) and the horizontal lines shows 
the false alarm probability at 1 and 2 a confidence limit. 

als from a linear ephemeris. Notice that the peaks close to the rota- 
tion period of the star have decreased. The main difference between 
this figure and the results for model 2 is due to the "spot features" 
during egress found for epochs 20, 21 and 147. Given that these 
"spot features" add power at the rotation period of the star, this 
suggests that these spots are real rather than systematic noise. 

However, we caution the reader to an over interpretation of 
these results because our model is very simplified, the fits are de- 
generate, and we cannot confirm that all the features present in the 
light curve are real spots. A complete spot model of WASP- 10 is 
beyond the scope of this paper, our objective being to estimate the 
mid-transit timing errors that can be caused by spots. 

In summary, stellar spots or correlated noise can significantly 
bias the derived parameters of the light curves especially if the spot 
occurs during ingress or egress. However, the errors due to spots 
can be accounted for if we include red noise due to shape variations 
as done for model 2. 

5.4 Effects of partial transits 

Several authors have noted that the uncertainti es in the TPs derived 
from partial trans it are underestimated (e.g iGibson et al.l 1200^ ; 
iBarros et al.l I2012T) . Given that the one transit that shows signifi- 
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Figure 12. Spot model transit residuals of WASP-10 after subtracting a linear ephemeris. The symbols are the same as Figure[6] RISE data are shown as filled 
stars, FTN data as squares while preciously published data as circles. 



cant transit timing variations is a partial transit (epoch 143), in this 
section we illustrate how the derived transit times are affected in 
partial transits. We assume here that the light curves are fitted si- 
multaneously (model 2), and hence the shape of the transits is set 
by the complete transits, i.e. we fixed a/R, and i in our test fits. 
To test the TTVs induced by partial transits, we used some of our 
complete light curves, from which we subtract model 2 estimated 
TO ( Table|6j. The light curves were then gradually cut at a time -Tp 
or Tp which corresponds to cutting the beginning or end of the ob- 
servation respectively. Resulting partial light curves were fitted to 
derived the induced TO variation. Red noise was estimated for each 
fake light curve and included in the errors. In Figure [T4l we show 
our test applied to the light curves epoch 0, 145, and 147 as exam- 
ples. For other light curves the derived TTVs are enclosed in the 
examples presented. We conclude that light curves that do not have 
red noise such as epoch show partial transit induced TTVs that 
are consistent with zero. However, for light curves with high red 
noise the induced TTVs can reach several minutes and the derived 
errors are indeed underestimated. The errors decrease with Tp due 
to red noise underestimation. For example, the TTVs derived from 
fake partials of epoch 147 are significant at 5 a even when we in- 
clude red noise. This is caused by the correlation between red noise 
and the normalisation parameters F ou t and T gra d in partial light 
curves. Notice that TTVs are induced even when Tp is larger than 
half of the transit duration, i.e. the transit needs to include some out 
of transit data. Therefore this could explain the 2.7 minutes TTV 
found for epoch 143, though we remind that this is also the only 
transit taken with 2.0m Rozhen telescope. 

If we discard all partial transits, the transit residuals of model 
2 are consistent with a linear ephemeris, with a \ 2 = 30 for 17 
degrees of freedom ( x? c d = 1-8) and a p- value = 0.025. Model 1 
results are still inconsistent with a linear ephemeris at the 5a level 
since transits show significant transit time variations due to spots 
that are not accounted for in this model as was shown in the previ- 
ous section. The one spot model timing residuals of the complete 
transits are also consistent with a linear ephemeris: \ 2 — 20 for 
17 degrees of freedom (x? c d = 1-2) and a p-value = 0.27. The 
increase of the p-value when we discard partial transits indicates 



i ° 

H -i 




Epoch 145 



/ Epoch 

iriiiiiiiii U¥ti±i 




Epoch 147 




-2-10 1 2 

cut (hours) 

Figure 14. Induced TTVs as a function of the time the light curves were 
cut, Tp. We show the errors accounting for red noise estimated for each fake 
partial light curve. The errors of epoch 145 were excluded for clarity. These 
errors are higher than for epoch 147, and hence, the induced TTVs are less 
significant. The vertical lines show the position of ingress and egress. 



that errors in transit times from partial transits are indeed underes- 
timated, as discussed above. 

In summary, it is possible to obtain a good estimation of the 
mid-transit times of partial transits if they do not have red noise. 
However, in the presence of red noise the errors of the mid-transit 
times of partial transits can be significantly underestimated due to 
our inability to recognise red noise in partial light curves . Unfortu- 
nately, complete transits are sometimes unattainable due to weather 
and technical issues. For example, all scheduled RISE observations 
were intended to cover the complete transit. We conclude that tran- 
sit times derived from partial transits have to be regarded with ex- 
treme care. 
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6 CONCLUSION 

We have presented eight new transit observations of WASP- 
10b with RISE mounted on the Liverpool telescope and Merope 
mounted on the Faulkes Telescope North. The new observations 
were combined with 22 previously published transit light curves to 
derive the transit times homogeneously. System parameters were 
re-derived assuming a circular orbit and assuming that the star is 
less active at maximum brightness. However, due to bright spots 
some stars, like the Sun, are actually brighter at the activity maxi- 
mum. No bright spot is seen in the transit cord of WASP- 10, there- 
fore, we assume that the observed stellar modulation is due to dark 
spots. 

We show that the significant transit timing v ariations for 
WASP- 10b reported by iMacieiewski et ail d201 lal fbh were due to 
the normalisation parameters being neglected in some of previous 
fits to the data. We find a time difference of 120 and 49 seconds for 
epochs 36 and 147, respectively when we account for a baseline 
function in the transit fit. It is crucial to account for the baseline 
function because it correlates strongly with the estimated transit 
times. 

We present two models to derive the transit times. Model 1 
where the transits are fitted individually and model 2 where the 
transits are fitted simultaneously. Comparison between the two 
models lead us to suspect that the transits times were affected by 
spot features or systematic noise. We show that spot features or 
systematic noise during ingress and egress have a significant effect 
on the measured transit times leading to spurious TTV detections. 
Since they also cause transit duration variations they are treated as 
red noise in model 2 which explains the differences found between 
the two models. Nevertheless, the transit times derived with model 
2 show higher dispersion than expected for a linear ephemeris. 

We tested if the transit times could be explained by a per- 
turber model. As expected the addition of a planetary companion in 
the system reduces the rms of the residuals. Ho wever, there is not 
more e vidence for the TTV signal claimed by Macieiews ki et al] 
neither in the periodogram nor in the n-body analysis 
of the TTVs. In fact, we found that other configurations with per- 
turber masses ranging from 0.01 to 5 M.j up and period between 1.5 
and 12 x the orbital period of WASP- 10b fit the data equivalently 
well. Therefore, we conclude that the transit times of WASP- 10b 
have insufficient accuracy and precision to validate the presence of 
any planetary companion. 

An alternative model for the dispersion seen in the transit 
times of WASP- 10b is stellar activity induced TTVs. We tested this 
hypothesis by including a spot feature in the transit model fit. Our 
simple model shows that transits 20, 21 and 147 are consistent with 
having a spot feature during egress that results in transit time vari- 
ations of -135, -78 and -126 seconds respectively. This explains 
the dispersion found for the transit times derived with model 2. 
We suggest that this spot induce transit variations are due to real 
spots rather than systematics because they add power to the residu- 
als at periods close to the orbital period of the host star. Moreover, 
other spots are seen during transit, in particular epochs 242, 243 
and 253 which were studied in detail by Ma cieiewski et ai]( l2011bl) 
that concluded that they were consistent with real spots. When we 
account for these spot features, the derived transit times are consis- 
tent with a linear ephemeris with the exception of epoch 143, which 
we show it could be due to this transit being partial. 

We also show that partial transits may produce fake TTVs up 
to 3 minutes, and hence TOs derived from partial transits should 
be regarded with extreme care. When we discard partial transits 



the dispersion of the transit times relative to a linear ephemeris de- 
creases significantly both for model 2 results and for the one spot 
model results. 



7 ACKNOWLEDGEMENTS 

We thank Damian Christian, John Johnson, Jason Dittmann, Tereza 
Krejcova, and Gracjan Maciejewski for kindly sharing their light 
curves on WASP- 10 and for interesting discussions. We thank the 
anonymous referee for contributing to the improvement of this 
manuscript. We thank Tom Marsh for the use of the ULTRACAM 
pipeline. SCCB thanks John Southworth and Rodrigo Diaz for in- 
teresting discussions. The RISE instrument mounted at the Liver- 
pool Telescope was designed and built with resources made avail- 
able from Queen's University Belfast, Liverpool John Moores Uni- 
versity and the University of Manchester. The Liverpool Telescope 
is operated on the island of La Palma by Liverpool John Moores 
University in the Spanish Observatorio del Roque de los Mucha- 
chos of the Instituto de Astrofisica de Canarias with financial sup- 
port from the UK Science and Technology Facilities Council. 



REFERENCES 

Adams E. R., Lopez-Morales M., Elliot J. L., Seager S., Osip D. J., 
2010, ApJ, 714, 13 

Agol E., Steffen J., Sari R., Clarkson W., 2005, MNRAS, 359, 567 

Alonso R., Aigrain S., Pont R, Mazeh T, CoRoT Exoplanet Sci- 
ence Team 2009, in IAU Symposium Vol. 253 of IAU Sympo- 
sium, Searching for the secondary eclipse of CoRoT-Exo-2b and 
its transit timing variations, pp 91-96 

Ballard S., Fabrycky D., Fressin R, Charbonneau D., Desert J.-M., 
Torres G, Marcy G, Burke C. J., Isaacson H, Henze C, Steffen 
J. H, Ciardi D. R., et al. 201 1, ApJ, 743, 200 

Barros S. C. C, Pollacco D. L., Gibson N. P., Howarth I. D., 
Keenan F. P., Simpson E. K„ Skillen I., Steele I. A., 201 1, MN- 
RAS, 416, 2593 

Barros S. C. C, Pollacco D. L., Gibson N. P., Keenan F. P., Skillen 

I., Steele I. A., 2012, MNRAS, 419, 1248 
Bean J. L., 2009, A&A, 506, 369 

Boisse I., Moutou C, Vidal-Madjar A., Bouchy R, Pont R, 
Hebrard G, Bonfils X., Croll B., Delfosse X., Desort M., 
Forveille T, Lagrange A.-M., et al. 2009, A&A, 495, 959 

Boue G, Oshagh M., Montalto M., Santos N. C, 2012, MNRAS, 
422, L57 

Christian D. J., Gibson N. P., Simpson E. K, Street R. A., Skillen 
I., Pollacco D., Collier Cameron A., Joshi Y. C, Keenan F. P., 
Stempels H. C, et al 2009, MNRAS, 392, 1585 

Czesla S., Huber K. R, Wolter U, Schroter S., Schmitt 
J. H. M. M., 2009, A&A, 505, 1277 

Demarque P., Woo J., Kim Y„ Yi S. K„ 2004, ApJS, 155, 667 

Dhillon V. S., Marsh T. R., Stevenson M. J., Atkinson D. C, Kerry 
P., Peacocke P. T, Vick A. J. A., Beard S. M., Ives D. J., Lun- 
ney D. W., McLay S. A., Tierney C. J., Kelly J., Littlefair S. P., 
Nicholson R., Pashley R., Harlaftis E. T, O'Brien K., 2007, MN- 
RAS, 378, 825 

Diaz R. R, Rojo P., Melita M., Hoyer S., Minniti D., Mauas 

P. J. D., Ruiz M. T, 2008, ApJL, 682, L49 
Dittmann J. A., Close L. M., Scuderi L. J., Morris M. D., 2010, 

ApJ, 717, 235 

Eastman J„ Siverd R., Gaudi B. S., 2010, PASP, 122, 935 



I TVs in WASP -10b induced by stellar activity? 17 



Fukui A., Narita N., Tristram P. J., Sumi T., Abe R, Itow Y., Sul- 
livan D. J., Bond I. A., Hirano T., Tamura M., Bennett D. P., 
Rurusawa K., et al. 201 1, PASJ, 63, 287 

Gelman A., Rubin D., 1992, Statistical Science, 7, 457 

Gibson N. P., Pollacco D., Simpson E. K., Barros S., Joshi Y. C, 
Todd I., Keenan R P., Skillen I., Benn C, Christian D., Hrudkova 
M., Steele I. A., 2009, ApJ, 700, 1078 

Gibson N. P., Pollacco D., Simpson E. K., Joshi Y. C., Todd I., 
Benn C, Christian D., Hrudkova M., Keenan F. P., Meaburn J., 
Skillen I., Steele I. A., 2008, A&A, 492, 603 

Gibson N. P., Pollacco D. L., Barros S., Benn C, Christian D., 
Hrudkova M., Joshi Y. C, Keenan F. P., Simpson E. K, Skillen 
I., Steele I. A., Todd I., 2010, MNRAS, 401, 1917 

Holman M. J., Fabrycky D. C, Ragozzine D., Ford E. B., Steffen 
J. H., Welsh W. F, Lissauer J. J., Latham D. W., Marcy G. W., et 
al. 2010, Science, 330,51 

Holman M. J., Murray N. W., 2005, Science, 307, 1288 

Howarthl. D., 201 1, MNRAS, 413, 1515 

Hoyer S., Rojo P., Lopez-Morales M., Diaz R. F, Chambers J., 

Minniti D„ 201 1, ApJ, 733, 53 
Husnoo N., Pont F, Mazeh T, Fabrycky D., Hebrard G., Bouchy 

F. , Shporer A. , 20 1 2, MNRAS , 422, 3 1 5 1 
Johnson J. A., Winn J. N., Cabrera N. E., Carter J. A., 2009, ApJL, 

692, LI 00 

Kipping D. M„ 2010, MNRAS, 407, 301 

Krejcova T, Budaj J., Krushevska V., 2010, Contributions of the 
Astronomical Observatory Skalnate Pleso, 40, 77 

Kundurthy P., Agol E., Becker A. C, Barnes R., Williams B., 
Mukadam A., 2011, ApJ, 731, 123 

Lissauer J. J., Fabrycky D. C, Ford E. B., Borucki W. J., Fressin 
F, Marcy G. W., Orosz J. A., Rowe J. F, Torres G., Welsh W. F, 
Batalha N. M., et al. 2011, Nature, 470, 53 

Maciejewski G., Dimitrov D., Neuhauser R., Niedzielski A., 
Raetz S., Ginski C, Adam C, Marka C, Moualla M., Mugrauer 
M., 2010, MNRAS, 407, 2625 

Maciejewski G., Dimitrov D., Neuhauser R., Tetzlaff N., 
Niedzielski A., Raetz S., Chen W. P., Walter F, Marka C, Baar 
S., Krejcova T, Budaj J., Krushevska V., Tachihara K., Taka- 
hashi H., Mugrauer M., 2011a, MNRAS, 411, 1204 

Maciejewski G., Raetz S., Nettelmann N., Seeliger M., Adam C, 
Nowak G., Neuhauser R., 201 lb, A&A, 535, A7 

Mandel K., Agol E., 2002, ApJL, 580, L171 

Miller-Ricci E., Rowe J. F, Sasselov D., Matthews J. M., 
Kuschnig R., Croll B., Guenther D. B., Moffat A. F. J., Rucinski 
S. M., Walker G. A. H., Weiss W. W., 2008, ApJ, 682, 593 

Nesvorny D., Kipping D. M., Buchhave L. A., Bakos G. A., Hart- 
man J., Schmitt A. R., 2012, Science, 336, 1 133 

Nesvorny D., Morbidelli A., 2008, ApJ, 688, 636 

Queloz D., Bouchy F, Moutou C, Hatzes A., Hebrard G, Alonso 
R., Auvergne M., Baglin A., Barbieri M., Barge P., Benz W., 
Borde P., et al. 2009, A&A, 506, 303 

Sanchis-Ojeda R., Winn J. N., Holman M. J., Carter J. A., Osip 
D. J., Fuentes C. I., 2011, ApJ, 733, 127 

Seager S., Mallen-Ornelas G, 2003, ApJ, 585, 1038 

Smith A. M. S., Hebb L., Collier Cameron A., Anderson D. R., 
Lister T. A., Hellier C, Pollacco D., Queloz D., Skillen I., West 
R. G., 2009, MNRAS, 398, 1827 

Steele I. A., Bates S. D., Gibson N., Keenan F, Meaburn J., Mot- 
tram C. J., Pollacco D., Todd I., 2008, in Society of Photo- 
Optical Instrumentation Engineers (SPIE) Conference Series 
Vol. 7014 of Edited by McLean, Ian S.; Casali, Mark, RISE: 
a fast-readout imager for exoplanet transit timing, pp 70146J- 



70 146 J 

Steffen J. H., Agol E., 2005, MNRAS, 364, L96 

Steffen J. H., Ragozzine D., Fabrycky D. C, Carter J. A., Ford 
E. B., Holman M. J., Rowe J. F, Welsh W. F, Borucki W. J., 
Boss A. P., Ciardi D. R., Quinn S. N., 2012, Proceedings of the 
National Academy of Science, 109, 7982 

Stetson P. B., 1987, PASP, 99, 191 

Winn J. N., Holman M. J., Torres G., McCullough P., Johns-Krull 
C, Latham D. W., Shporer A., Mazeh T, Garcia-Melendo E., 
Foote C, Esquerdo G., Everett M„ 2008, ApJ, 683, 1076 



